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ABSTRACT 

The Cygnus Loop has been the focus of substantial debate concerning the contri¬ 
bution of charge exchange (CX) to supernova remnant (SNR) X-ray emission. We 
take advantage of a distinct feature of CX, enhanced Ka forbidden line emission, and 
employ the energy centroid of the OVII Ka triplet as a diagnostic. Based on X-ray 
spectra extracted from an extensive set of Suzaku observations, we measure the energy 
centroid shifts of the triplet on and off the shock rim of the remnant. We find that 
enhanced forbidden to resonance line emission exists throughout much of the rim and 
this enhancement azimuthally correlates with non-radiative Ha filaments, a tracer of 
strong neutral-plasma interaction in the optical. We also show that alternative mech¬ 
anisms cannot explain the enhancement observed. These results demonstrate the need 
to model the CX contribution to the X-ray emission of SNRs, particularly for shocks 
propagating in a partially neutral medium. Such modelling may be critically impor¬ 
tant to the correct measurements of the ionization, thermal, and chemical properties 
of SNRs. 

Key words: atomic processes; ISM: supernova remnants; ISM: abundances; ISM: 
individual: Cygnus Loop; X-rays: ISM 


1 INTRODUCTION 


X-ray astrophysics is layered with a host of interesting phe¬ 
nomena. As spectral capabilities continue to grow within 
the X-ray regime, our ability to discern the importance of 
these varied phenomena has led to many interesting sur¬ 
prises. One X-ray emission mechanism whose importance 
we are increasingly becoming aware of is charge exchange 
or charge transfer (CX; e.g., Lallement 2004). While CX re¬ 
actions have been studied in an astrophysical context for 
decades, their importance for X-ray astro physics wasn’t re¬ 
alized until c omet observations in 1996 (iLisse et al.l 1 19961 : 
ICravensl fl^. Since this cometary epiphany, CX has been 
proposed and studied in a wide range o f X-ray sources, from 
planetary atmospheres (iDennerll l2002l : iDennerl et al.l l2006l : 
[Kharche nko et ^ |2008 |) throug h individual star-forming 
complexes (|Ttawrisley_etLS- 2011) and s tar-forming galaxies 
(iTsum et al.ll200';1:lLiu et al.ll201ll.l2012l) to the cores of clus¬ 
ters I Fabian et al ][2on|). More recently, the first evidence for 
X-ray emit ting CX in a superno va remnant (SNR) has been 
suggested llKatsuda et al.ll2011^ . 

CX is a special case of inelastic collision in which one or 
more electrons is exchanged between an ion and typically a 
neutral atom, most likely hydrogen. As the colliding parti¬ 
cles approach each other, their potentials overlap, allowing 
the photonless exchange and subsequent cascade of the cap¬ 
tured electron. This is an extremely efficient process with 
a relatively large reaction cross section (about 5 orders of 


* E-mail: srrobert@astro.umass.edu 


magnitude greater than for electron collisional excitation). 
If the capturing ion is sufficiently ionized (e.g. OVIII) then 
the cascading is expected to result in X-ray emission. Thus, 
in principle. X-ray emitting CX can be an important pro¬ 
cess at any astrophysical site where heavily ionized species 
interact with neutrals, making the thin post-shock region of 
SNRs a very promising site for CX. 

While CX has been used to explain the broad compo- 
nent of Ha emission in many SNRs for more than 30 years 
llChevalier et al.lll980l : [Chavamian et al.ll200ll) . X -ray emis¬ 
sion produced by CX has not yet been directly apparent. 
This is at least partly due to the lack of a non-dispersive 
high-resolution spectrometer. We are instead confined to X- 
ray CCD data for our studies of diffuse emission. At the 
spectral resolution of such data, the pure line emission of 
CX can mimic the continuum + line emission of a thermal 
plasma, making it difficult to disentangle the relative contri¬ 
butions of the different emission mechanisms. Thus, without 
accounting for the contribution from CX, the modelling of 
the properties of the shocked plasma (e.g. chemical abun¬ 
dances) can potentially be misleading. 

To delve into the importance of CX for SNR emission 
we need to be able to spatially and spectrally resolve the 
nuances of the plasma evolution in nearby remnants, such 
as the Cygnus Loop SNR. This remnant has an angular di¬ 
ameter of « 3° and suffers relatively little absorption. Such 
a close SNR provides a unique laboratory for the study of 
shock physics and non-equilibrium conditions. These spatial 
and spectral studies have been done to some extent in the 
X-ray regime with Suzaku and XMM-Newton, for which ob- 
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c 

N 

O 

Ne(=Mg) 
Si(=S) 
Fe (=Ni) 


0.164 (±0.016) 
0.130 (±0.020) 
0.105 (±0.009) 
0.230 (± 0.016) 
0.243 (±0.043) 
0.125 (±0.011) 


Table 1. Typical best fit metal abundances (fraction of solar) for 
the immediate post shock region, taken from the NE4 observation 
(see Section [^for more details). 


servations cover nearly the entire remnant. However, these 
studies have seemingly created more questions than answers. 

With Suzaku, spectral fitting of the immediate post¬ 
shock region has led to a quagmire of bewilderment where 
derived metal abundances appea r to fall into two gen¬ 
eral categories (jUchida et al.ll2f)^ ). In some of the rim re¬ 
gions, abundances have been cited as 0.5-1 times the so- 
lar value, which is consistent with the surrou nding ISM 
l|CartIedge et al.l l200fil : iNieva fc PrzvbiU^ 120121 ) . Through¬ 
out the rest of the rim, abundances are found to be about 
0.2 times the solar value or less, well outside the dispersion 
of local ISM measurements (e.g. Table [T]and Figure [T]for 
the results of a typical spectral fitting in the immediate post¬ 
shock region). This peculiar situation has led to the consid¬ 
eration of phenomena that are either generally thought to 
be unimportant or more exotic and poorly understood, with 
no clear answer emerging. 

Missing physics in the modelling of moderate resolu¬ 
tion spectra is probably the most natural and intuitive 
explanation f or the regions of apparent low abundance. 
iMivata et al.l (l2008l ) suggested that resonance-line scatter¬ 
ing might account for the regions of relatively suppressed 
abundances. However, the authors found that even with as¬ 
sumptions that are likely to amplify the effects of scatter¬ 
ing, it could not account for the level of abundance suppres¬ 
sion found, particularly that of o xyge n. Reasoning along the 
same lines. [Katsuda et al.l ll2008l) an d iTsunemi et al.l (l2009l ) 
posited that the continuum could be higher in these regions 
due to a nonthermal contribution, thus creating the appear¬ 
ance of lower metal abundances. Although this did improve 
the quality of the fits and increased the abundances, the pre¬ 
dicted radio fluxes are orders of magnitude higher than the 
observed values. 

With strong line emissio n, CX could significant ly affect 
the abundance measurement. iKatsuda et al.l |20l3) noticed 
when excluding energies around 0.7 keV, the abundances 
typically drop to those of the low abundance regions. They 
posited that the high abundances may be artificial, proba¬ 
bly due to CX contributions at ~ 0.7 keV from higher order 
transitions of O vii (note the excess emission around this 
energy in Figure [!}. Also, this enhancement is found to 
anti-correlate with non-thermal radio emission. This led the 
authors to argue that the cosmic-ray precursor more effec¬ 
tively ionizes the pre-shock ISM, thus eliminating the CX 
process. However, if this interpretation were to prove valid, 
it would also necessitate an explanation for the overall low 
abundances throughout the rim. Furthermore, the rim re¬ 
gions where the enhancement occurs do not correlate with 
non-radiative Ha emission (an indicator of neutral gas in 
the post-shock plasma), as would be expected for CX. 

Here we propose to judge the importance of CX based 
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Figure 1. Typical spectrum of the immediate post-shock region, 
taken from the NE4 observation (see Section [2] for more details). 
Spectrum is fit with an absorbed single temperature NEI plasma 
with best fit abundance parameters listed in Table [T] 


on a less model-dependent diagnostic. In this work we will 
exclusively use the O vii Ka triplet (-^ 0.57 keV), which 
is often the strongest line in the spectra. The triplet con¬ 
sists of the resonance (574.0 eV), inter-recombination (568.6 
eV) and forbidden (561.1 eV) transitions. Effects of CX 
become clear when the so-called G ratio [= (/ ± i)/r, 
forbidden±inter-combin ation to resonance t ransitions] of 
the triplet is analyzed (iLiu et al.l l201ll . |2012|) . Due to the 
cascading and high number of quantum states leading to 
the forbidden line, the ratio is calculated to be enhanced 
relative to thermal emission (due to collisional excitation). 
Even if unable to resolve the individual components of the 
triplet, one may adopt its central energy as a proxy for the 
G ratio, which can typically be determined to an accuracy of 
more than 20 times better than the intrinsic spectral resolu¬ 
tion when the signal to no ise ratio is sufficiently high, which 
is typically the case here (IWargelin et al.ll2008l) . Therefore, 
the centroid analysis can be useful, even with a spectrum of 
moderate resolution. For the O vii Kq triplet, the centroid 
is expected to be shifted to lower ene rgies by 4.3(±0 .5) eV 
for GX emission compared to thermal (ISnowderJl2009l ) . This 
enhancement in the forbidden line is illustrated in Figure [T] 
by the large residuals on the low energy side of the O vii Ka 
line. 

Our observations and data analysis are described in § 
[2] We explore how the central energy of the triplet changes 
azimuthally and with distance from the shock front in § |31 
The spatial distribution of the energy shift is then compared 
directly to a tracer of strong plasma-neutral interaction in 
§ m Lastly, we compare these results with previous studies 
of the Cygnus Loop SNR and discuss the implications and 
alternative scenarios in § (5] 


2 OBSERVATIONS AND DATA REDUCTION 

Suzaku, the ideal telescope for this study, was launched 
on July 10, 2005 carrying 4 CGDs with independent X- 
ray telescopes sensitive to photons in the 0.2-12 keV range 
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Figure 2. Cygnus Loop SNR regions observed with Suzaku (17.8^ 
X 17.8') overlaid on the logarithmically scaled Rosat HRI im¬ 
age. Cyan observation regions do not show an energy shift of the 
O VII Ka line from the rim to the inner region while black ones 
do (see text for details). 


ijKunieda et al.l [20061 ). While three of the CCDs are front- 
illuminated, one (XISl) is back-illuminated and is more sen¬ 
sitive to soft X-rays. Thus, the XISl spectrometer (with a 
spectral resolution of ~ 60 eV at the energy of the OVII Ka 
triplet) is used for the present study. Observations of the 
Cygnus Loop with Suzaku encompass nearly the entire rim 
of the SNR, excluding the region of the southern blowout 
and a small area in the eastern face of the shock front (Ta¬ 
ble [21 Figure (21). 

The observations are reduced according to the standard 
procedures described in the Suzaku ABC guide, using HEA- 
soft version 6.11.1 and CALDB version 2011-11-09. The data 
is reprocessed using the ftool aepipeline. When available, 
data from the 5x5 and 3x3 editing modes are combined. Two 
selections are placed on the data, cut-off rigidity (COR )>8 
and elevation from Earth (ELV)>10°. The COR selection is 
used to limit the contamination from the non-X-ray back¬ 
ground (discussed below). The elevation selection is made 
to keep the Oxyge n 0.54 keV flourescent line to a minimum 
llSmith et al.ll200% . which would be blended with the OVII 
Kq triplet, potentially leading to problematic results. The 
channels of all spectra are adaptively grouped to have a min¬ 
imum (background-subtracted) S/N threshold of 3 per bin 
before fitting. 

Since our study is based on the differential line centroid 
of the OVII Kq triplet on and off the SNR shock, it is im¬ 
portant to check whether or not any systematic instrumental 
gain shift may affect our measurements. Investigation of the 
CCD shows no apparent systematic shift along the CCD. 


Thus, we can compare the relative line centroids between 
regions within an observation without having to consider 
the actual gain error of the instrument. 


3 SPECTRAL ANALYSIS 

Our spectral analysis is conducted with the XSPEC software 
package (version 12.7.1). In order to characterize the shock 
region of the Cygnus Loop SNR, spectra are extracted from 
2’ wide cuts along the shock front. Following the findings 
of llKatsuda et al.l 201ll ~l. in which the authors determined 
any enhancements at ~ 0.7 keV occurred within 4’ of the 
shock front, spectra were extracted from the innermost re¬ 
gions of the same observations up to 4’ from the shock front 
in order to characterize the inner regions of the SNR. Such 
a spatial dist ribution is also c onsistent with previous theo¬ 
retical work ([Lallementl l2004h . These two spectra are then 
compared directly to each other. 

The background can arise from three general sources: 
extragalactic, particle, and locally (SXB). The particle back¬ 
ground is modelled using night Earth observations with the 
ftool xisnxbgen and is subtracted from the source spectra in 
XSPEC before spectral fitting. The SXB contribution due to 
the Local Bubble and solar wind CX at Earth’s outer atmo¬ 
sphere and the heliosphere, is modelled as an unabsorbed 
single temperature thermal plasma, which has been found 
to give adequate fits de spite the contested importance of 
CX to the SXB emission ([Anderson et al.ll2010l : [Smith et al.l 
l2014h . The temperature of this local plasma component 
is hxed to the fiduci al value previously found, 0.11 keV 
([Anderson et al.ll2010[ l . This local emission is dwarfed by the 
SNR emission and contributes only about 1% of the emission 
in soft X-rays. The extragalactic component, which domi¬ 
nates the emission at high energies after particle background 
remova l, is modelled as a p ower law with a photon index of 
r=1.4 ([Yoshino et al.[l2009[f together with a f oreground ab¬ 
sorpt ion column density of 3 x 10^° cm“^ ([Kosuei et al.l 

l2010[f . 

In order to characterize the center of the O vii Kq line, 
it must be modelled as a separate gaussian. First, we mod¬ 
elled the source emission as a non-equilibrium ionization 
(NEI) plasma with floating abundances of each element indi¬ 
vidually and foreground absorption column density frozen at 
3 X 10^° cm~^. After this initial fit is obtained, the oxygen 
abundance is set to 0 and all other parameters are frozen. 
The Oxygen emission is then modelled as a set of lines: OVII 
Kq, OVIII LyQ, OVII K/3, OVII K 7 , OVII Kd, OVII Ke, 
OVIII LyP, OVIII Ly 7 , OVIII Lyd, and OVIII Lye (all line 
energies hxed to theoretical val ues except for the O v ii Kq 
and O VIII LyQ lines) as done in iKatsuda et all ([201 il l. The 
lines are not expected to be resolved and are set to have a 
width of 0 keV. 

The difference in the OVII Kq centroid energy between 
the inner and rim spectra is presented in Table [3] for each 
observation. While several of the observations are consis¬ 
tent with no shift, most of the observations do exhibit an 
apparent shift as expected for CX emission. Comparing the 
observations that have an apparent centroid shift (i.e., > 3 
eV) to Figure (2) we see that nearly the entire rim of the 
SNR sees a shift to higher energies from the rim inward with 
the exception of two regions, one in the northwest rim and 
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Table 2. Suzaku XISl Observations of the Cygnus Loop SNR 


Obs ID 


1 

b 

Start 

End 

to 

te 

500020010 

(NEl) 

75.61 

-8.84 

2005-11-23 17:39:01 

2005-11-24 04:55:24 

20.34 

14.36 

500021010 

(NE2) 

75.66 

-8.56 

2005-11-24 04:56:05 

2005-11-24 16:14:24 

21.41 

20.88 

500022010 

(NE3) 

75.72 

-8.28 

2005-11-29 17:47:47 

2005-11-30 05:39:09 

21.13 

17.34 

500023010 

(NE4) 

75.72 

-7.99 

2005-11-30 05:41:02 

2005-11-30 18:23:14 

25.00 

19.02 

501020010 

(PIO) 

75.13 

-7.94 

2007-11-13 06:01:00 

2007-11-13 14:30:15 

16.80 

10.19 

501035010 

(P18) 

72.85 

-8.69 

2006-12-18 01:16:09 

2006-12-18 08:10:21 

12.01 

8.71 

501036010 

(P19) 

73.01 

-8.29 

2006-12-18 08:11:18 

2006-12-18 19:47:24 

18.60 

17.92 

503057010 

(P21) 

75.60 

-7.76 

2008-06-02 22:33:53 

2008-06-03 07:08:24 

16.17 

10.24 

503058010 

(P22) 

75.39 

-7.53 

2008-06-03 07:09:17 

2008-06-03 18:03:14 

19.29 

16.49 

503059010 

(P23) 

75.17 

-7.32 

2008-06-03 18:04:23 

2008-06-04 03:32:14 

19.46 

15.46 

503060010 

(P24) 

74.93 

-7.14 

2008-06-04 03:33:03 

2008-06-04 15:03:09 

18.50 

13.50 

503061010 

(P25) 

74.69 

-7.04 

2008-06-04 15:04:06 

2008-06-05 03:34:24 

26.00 

21.16 

503062010 

(P26) 

74.46 

-9.71 

2008-05-13 02:53:19 

2008-05-13 13:08:19 

16.92 

14.62 

503063010 

(P27) 

74.06 

-9.70 

2008-05-13 13:09:08 

2008-05-14 01:11:14 

22.78 

16.94 

503064010 

(P28) 

73.77 

-9.54 

2008-05-14 01:12:11 

2008-05-14 12:48:14 

18.17 

15.72 

504005010 

(Riml) 

74.35 

-7.07 

2009-11-17 22:40:59 

2009-11-18 22:46:24 

40.75 

27.30 

504006010 

(Rim2) 

74.01 

-7.10 

2009-11-18 22:47:37 

2009-11-19 11:37:24 

26.31 

14.89 

504007010 

(Rim3) 

73.12 

-7.21 

2009-11-19 11:38:17 

2009-11-20 02:53:24 

21.56 

17.63 

504008010 

(Rim4) 

73.57 

-7.49 

2009-11-20 02:54:21 

2009-11-20 08:35:24 

12.10 

5.63 

504009010 

(Rim5) 

73.33 

-7.73 

2009-11-20 08:40:05 

2009-11-20 19:05:17 

15.85 

12.05 

504010010 

(Rim6) 

75.49 

-9.15 

2009-11-20 19:10:13 

2009-11-21 04:11:24 

14.33 

10.44 


Table 3. OVII Ka centroid line difference between rim and inner 
spectra {Eisner — Erim)- Errors are at the Icr confidence level. 


Obs ID 

AE (eV) 

501020010 (PIO) 

6.0 (± 

1.41) 

501035010 (P18) 

4.9 (± 

1.35) 

501036010 (P19) 

3.1 (± 

0.67) 

503057010 (P21) 

5.1 (± 

0.86) 

503058010 (P22) 

3.8 (± 

0.58) 

503059010 (P23) 

3.7 (± 

0.82) 

503060010 (P24) 

3.9 (± 

1.05) 

503061010 (P25) 

4.0 (± 

1.26) 

503062010 (P26) 

0.9 (± 

0.92) 

503063010 (P27) 

-1.0 (± 

0.95) 

503064010 (P28) 

-2.0 (± 

0.79) 

504005010 (Riml) 

-0.9 (± 

1.39) 

504006010 (Rim2) 

-0.9 (± 

1.04) 

504007010 (Rim3) 

0.0 (± 

0.93) 

504008010 (Rim4) 

5.0 (± 

1.41) 

504009010 (Rim5) 

8.0 (± 

1.83) 

504010010 (Rim6) 

5.0 (± 

1.12) 

500020010 (NEl) 

5.0 (± 

0.77) 

500021010 (NE2) 

3.0 (± 

0.58) 

500022010 (NE3) 

6.0 (± 

0.35) 

500023010 (NE4) 

6.0 (± 

0.34) 


another in the southeast rim. These two regions also have 
the lowest downstream X-ray flux, with the exception of the 
blowout, due to the relatively lower density of the ISM in 
these regions. The lack of a centroid shift in these regions is 
probably a natural consequence of a shock propagating into 
a low-density medium, which also presumably has a rela¬ 
tively small neutral hydrogen fraction and hence little CX 
emission is expected. 


4 COMPARISON WITH NON-RADIATIVE Ha 

While radiative Ha emission is associated with radiative 
cooling, non-radiative Ha is associated with relatively faster 
shocks (vs ^ 200km/s) that heats the largely neutral up¬ 
stream gas to a temperature that emits X-ray emission. 
Such shocks are characterized by the lack of emission from 
forbidden lines of lowly-ionized metals. This non-radiative 
Ha emission arises primarily from hot protons after electron 
capture via CX with neutral atoms and from collisional exci¬ 
tation of cold neutral hydrog en in the immediate post-shock 
region (iLevenson et al.[|l99^ . With spectroscopy these two 
contributions can be separated due to velocity distribution 
differences of the two sources with CX producing a broad 
line and collisional excitation producing a narrow line. How¬ 
ever, both components are related to the interaction of neu¬ 
trals with the hot plasma behind the shock front. Therefore, 
the Ka emission of OVII via CX should be correlated with 
regions of enhanced Ha emission relative to low ionization 
metal forbidden lines. 

In order to compare with non-radiative Ha emission 
of the Cygnus Lo op SNR, we use the g enerously provided 
data presented in ILevenson et ahl (Il998ll . However, the Ha 
emission can arise from both radiative (cooling) and non- 
radiative regions. To extract the structure and emission that 
is predominantly non-radiative, we use the ratio of two nar¬ 
row band images, Ha to [SH] (AA6717,6731) (IR; Figure |3]). 
Notice that the regions dominated by Ha are confined pre¬ 
dominantly to the rim of the SNR. Since this emission is 
relatively weak and only occurs in a narrow region behind 
the shock, similar to the CX X-ray emission, it only becomes 
apparent where the path length through the emission region 
is long. Visual inspection indicates that there is a clear cor¬ 
respondence with the regions of the centroid shift and we 
see that strong non-radiative filaments are lacking in the 
southeast and northwest regions of the SNR. 

To quantify this correlation, we extract the Suzaku ob- 
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Figure 3. Logarithmically scaled Ho/[S II], 5R, ratio map 
smoothed with a 3 pixel kernel. Green is equivalent to a 5R = 1, 
with red being the regions of large SR (relatively high Ho) and 
blue being the regions of small SR. The red contour shows the 
minimum ratio value to be considered non-radiative Ho, which 
we see are predominantly confined to the rim of the SNR. Data 
courtesy Dr. Nancy Levenson (Levenson et al. 1998). 



Figure 4. OVII Ko centroid shift between inner and rim regions 
(Einner — Ej-im) as a function of integrated Ho flux density (nor¬ 
malized to 1) in observation region. All data points of relative 
integrated H-alpha flux below 0.2 are in actuality at 0. These 
points are staggered to allow one to see the different points. We 
find that the presence of non-radiative Ha filaments correlates 
with a centroid shift with 97% certainty. Reported error bars are 
1 O'. 

was left out of the analysis due to a saturated point source, 
52 Cygni, in the FoV of the unsmoothed optical image. 


servation regions individually from the unsmoothed global 
Ha to [S II] ratio map. Non-radiative H-alpha emission is 
then found by computing dendrograms of the image with 
minimum threshold K„iin=l.l and a minimum number of 
pixels of 200, which was chosen by eye to only select true 
regions of filaments. The algorithm selects out continuous 
structure starting from the pixel maxima with the given 
threshold criteria. This selection criteria is highlighted by 
a red contour in Figure [3] 

All structures found in the extracted regions are then 
combined. These structures found via K are then mapped 
to Ha images and the integrated flux from the structures 
is calculated. A comparison of this integrated flux with the 
centroid shift is shown in Figure |4] We see that observations 
without a centroid shift are confined to lower Ha hlament 
flux, with little filamentary structure. Additionally, these X- 
ray centroid shift measurements appear to separate into an 
apparent dichotomy, those that have no shift and those that 
are forbidden line dominant and saturated. While the deter¬ 
mination of whether this saturation is physically reasonable 
is beyond the scope of this paper, we can ask the simple ques¬ 
tion of whether the presence of non-radiative Ha filaments 
correlates with an OVII Ka centroid shift, as we might ex¬ 
pect if the soft X-ray emission is significantly contaminated 
by CX emission. 

Binning the data based on whether there is an apparent 
centroid shift and on the presence of non-radiative Ha fila¬ 
ments, we use Fisher’s exact test to get a statistical handle 
on whether these two phenomena are correlated. We cal¬ 
culate a p-value of 0.03, indicating significant likelihood of 
association between an OVII Ka centroid shift and the pres¬ 
ence non-radiative Ha emission. Note, the Rim5 observation 


5 DISCUSSION 

5.1 Comparison with Previous Work 


Radially both this analysis and that of iKatsuda et ah ^2011 ) 
agree with the theoretical calcuations done by iLallementl 
ll2004l) . at least qualitatively. In that work, the author shows 
that CX effects are generally negligible for an SNR in the 
Sedov phase, except in the narrow rim region where line 
of sight effects amplify the fraction of emission from CX. 
In Figure 3 of that work, it is shown that CX can con¬ 
tribute over 50% of the X-ray emission immediately behind 
the shock and approximately 30% in the re gion contained 
withi n the shock radius(Rs) and r=0.98Rs i Katsuda et al.l 
l201ll) . thus detailing the confinement of CX to the rim. How¬ 
ever, azimuthally we find that the regions with OVII Ka 
line centroid shifts do not match those wit h the emission 
enhan cement at 0.7 keV, as studied by IKatsuda et al.l 
1I2OIII) . In fact, all of the regions that do not show an excess 
at ~ 0.7 keV exhibit a centroid shift. 

In overall flux, ther e is a great deal of ten sion between 
the CX interpretati on ofiKatsuda et al.l (I2OIII) and the the¬ 
oretic al results of IWise fc Sarazin 1 19891) . Katsuda et al.l 
1I2OIII) estimate the total CX contribution to the overall 
SNR emission to be 1- 3 orders of magnitude higher than 
the prediction made by IWise fc SarazinI lll989l) ; the authors 
note, however, that the theoretical calculations were done 
with different shock velocities than those of the Cygnus Loop 
S NR. There is also slight tension between our result and that 
of lWise fc SarazinI (Il989l ) . In that work they predict the CX 
contribution to the OVII Ka line to be at most about 1% 
of the line flux over the entire remnant. If we assume that 
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all of the line emission in the rim regions are due to CX 
and that the bulk of the emission of the SNR occurs in the 
outskirts of the SNR (within our observations shown), we 
roughly estimate the CX contribution to the O vii Ka line 
to be « 10%, as an upper limit. 

Beyond the reasonable agreement with theoretical work, 
our CX scenario is further confirmed by the correlation with 
non-radiative Hq filaments. Also, this analysis is far less 
model dependent. It seems natural that this contamination 
of the X-ray spectra by CX may lead to the underestimate of 
the metal abundances when they are modelled with an over¬ 
simplified plasma model. However, the lack of spatial corre¬ 
lation between the OVII Kq shifts and 0.7 keV enhancement 
(which traces the abundance inhomogeneities) suggests that 
CX is not the only phenomenon miss ing in the modelling . 
This is the same conclusion reached bv ICumbee et all (|2014l 'l 
where the authors attempted to do some modelling of CX 
emission in the remnant. 


5.2 Resonance-line scattering 

There are a couple other mechanisms that could serve to in¬ 
crease the G ratio to the levels expected by CX. Resonance¬ 
line scattering is one such mechanism. As seen in this work, 
the centroid energy shift is confined to the thin post-shock 
region. This is also expected from resonance-line scattering 
as the optical depth may be greatly enhanced along tangen¬ 
tial lines of sight through the immediate post-shock region 
where density is the greatest. Resonance line scattering has 
been shown to be potentially important, causing the line in¬ 
tensity of the oxygen Ka triplet being underestimated by up 
to 40% in regions th at show apparently low metal abundance 
l|Mivata et al.ll20(^ 'l. A reduction in the triplet intensity of 
this magnitude due to resonance-line scattering would result 
in an increase of the G-ratio to five assuming an initial G- 
ratio of one, as is typical of (nearly) thermal emission, which 
would be more than enough to shift the line energies to the 
extent that we have measured. 

We examine the importance of the resonance scattering 
by confronting its expected and observed spatial and flux 
dependences. In many of the regions that exhibit a centroid 
shift, the emission measure is larger in the inner spectrum 
than in the rim, while also having a temperature and ioniza¬ 
tion timescale that favors OVII as the dominant ionization 
state. Consider a plasma of electron temperature equal to 
0.21 keV, as is typical of our post-shock fits. The main con¬ 
stituent throughout a wide range of ionization timescales 
would be OVII, from 1.x 10® s/cm® to 4.x 10®® s/cm®. Once 
the plasma is equilibrated, the OVII fraction is still approx¬ 
imately 0.250 Combining the mild change in OVII fraction 
with the larger interior emission measure, we would expect 
resonance line scattering to be more important interior to 
the rim. 

For example, in the NEl region the best fit emission 
measures are 5.62x 10®®^ cm“® in the rim {EMr) and l.OGx 
10®® cm“® in the inner region (EMi). Assuming the regions 


^ These calculations were done using AtomDB 3.0, which is cur¬ 
rently in a beta testing phase. Its atomic physics is updated to 
include the effects of inner-shell ionization up to Nickel. 


have the same intrinsic oxygen abudance, the ratio of their 
optical depths can be written as: 

fr ^ /oVII.rtte.r At- _ fo'Vll,rE]\d^ Lj- . . 

Ti ~ fovu,ine,iLi /ovii.iEM®/®!.®^® 

where fovii is the ionic fraction of OVII, Ue the electron 
density, and L is the path length through the plasma. As¬ 
suming a spherical remnant and using the path lengths at 
the inner sides of the extraction regions (~ 2’ and 6’) re¬ 
sults in Lr/Li ~ 0.53. Using a conservative fovii,T-/fovii,i 
~ 2, equivalent to fully equilibrated in the inner spectrum 
and typical fractions of OVII for our post-shock ionization 
timescale fits, results in a strict upper limit of Tr/rt ~ 1 .1. 

There are a couple reasons why this opacity ratio is 
likely to be an overestimate. First, if one fixes the metal 
abundances in the rim to the local ISM values, the best 
fit emission measure is reduced by factor of 10, effectively 
reducing the calculated optical depth by a factor of 3. Sec¬ 
ond, adding any missing line emission physics (e.g. GX) will 
likely serve to suppress the emission measure further, par¬ 
ticularly in the rim since we expect this region to be further 
from equilibrium and thus more likely to suffer from missing 
physics. Taken in conjunction, we expect the optical depth 
in the rim to actually be quite a bit lower than interior. 

Even in the most favorable case, the optical depths in the 
rim and inner regions are comparable and we do not expect 
the differential resonance scattering between regions could 
result in our measured line centroid shifts. 

5.3 Inner-shell ionization 

Another process, inner shell ionization, is characteristic of an 
NEI plasma and has been show n to be import ant in the spec¬ 
tral modelling of young SNRs (lKosenkdl2006li . Through this 
mechanism, the G-ratio can be increased when the plasma is 
underionized and OVI is the main constituent. When an OVI 
ion undergoes inner shell ionization an electron is ejected 
from the ground state, leaving one electron in the Is ®So 
state and the other in the ls2s ®So or ls2s ®Si state. This 
would typically lead to the forbidden transition and hence 
the G-ratio will be enhanced. Figure [5] demonstrates the 
potential effects of this process. 

In this figure we calculate the cumulative G-ratio for 
a slit aperture of increasing width from the shock front to¬ 
ward the center of the SNR for a spherical shock of the size 
of the Cygnus Loop SNR. The calculation is done spanning 
the range of surrounding ISM densities typicaUy cited us¬ 
ing updated atomic data provided by Borkowsko to include 
inner-shell ionization and assuming infinite spatial resolu¬ 
tion. Note that the calculation also includes contributions 
from lithium-like satellite lines («560eV), which are blended 
with the forbidden line. 

Lower ISM densities push the curve to the right by low¬ 
ering the ionization timescale at fixed distance from the 
shock front. The G-ratio quickly peaks at « 6.7, but de¬ 
creases sharply beyond the peak and asymptotically ap¬ 
proaches the collisional equilibrium value within several arc- 
seconds, even for the favorable conditions considered here. 

The state in which this effect is important occurs very early 

^ http://space.mit.edu/home/dd/Borkowski/APEC_nei_README.txt 
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Figure 5. The cumulative G-ratio plotted as a function of the 
angular width extending from the shock front when inner shell 
ionization of Li-like oxygen is included. Different lines correspond 
to different pre-shock gas densities, as specified in the insert, but 
assume the same characteristic shock velocity of 330 km/s and 
distance (540 pc) of the Cygnus Loop SNR. This shock velocity 
is calculated from the typical electron temperature found in the 
immediate post-shock region (0.21 keV). 


on, at net « 10® s/cm®, making this effect negligible for the 
angular resolution of the study presented here and unable 
to explain the O vii Ka centroid shift observed. 


6 CONCLUSIONS AND FUTURE PROSPECTS 

Plasma interfaces are regions rich in interesting physics. We 
are just now straining to visualize their intricacies as we push 
the limits of our observational ability. For instance, without 
being able to identify spectroscopic signatures of CX and 
inner-shell ionization, such as Ka G-ratios, the overall con¬ 
tributing effects have been elusive for X-ray astrophysics and 
it becomes difficult to properly infer the plasma properties, 
such as temperature and metal abundances, particularly for 
the CX case since it has yet to be well modelled. Ideally, one 
wants to resolve the Ka triplets and hence to use its full di¬ 
agnostic power. Nevertheless, we can still utilize unresolved 
Ka triplets in cases of either moderate resolution spectra or 
low S/N high-resolution spectra to probe the physical con¬ 
ditions of plasma via their centroid, as demonstrated in the 
presented analysis of the Cygnus Loop SNR. The conclu¬ 
sions of this analysis are as follows: 

• The centroid shift between the on and off-shock OVII 
Ka triplet can be naturally explained by CX in terms of its 
magnitude (« 5 eV) and radial distribution. 

• The azimuthal correlation of the shift with non- 
radiative Ha emission reinforces this interpretation. 

• CX appears to be the only viable mechanism to explain 
the centroid shift. 

• CX as traced by the centroid shift does not spatially 
correlate with, and therefore can’t be the sole source of, the 
abundance inhomogeneities. 


These conclusions seem to be firm, demonstrating the po¬ 
tential of this novel approach. Even as we move into the 
future with the next generation of X-ray telescopes, all with 
high-resolution spectroscopic capabilities, approaches simi¬ 
lar to the one presented in this paper will remain useful. 
Astro-H, for example, will have an X-ray calorimeter on¬ 
board, with an expected angular resolution better than 1.3 
arcmin and expected spectral resolution of ~ 7 eV. While 
such an instrument will undoubtedly be revolutionary for 
our understanding of diffuse X-ray emission, even with this 
resolution the triplet cannot be fully resolved, making cen¬ 
troid energy diagnostics vital to the future of X-ray spec¬ 
troscopy. Unfortunately, the small FoV of Astro-H (similar 
to the angular resolution), will prevent the specific approach 
presented here from being overly useful. However, one may 
reasonably use consecutive observations with little fear of 
a change in instrumental gain. Furthermore, CX is not ex¬ 
pected to influence the energy of the OVHI Ka line. This 
allows us to leverage this line as an instrumental gain ref¬ 
erence and probe any OVII forbidden line enhancement for 
the case where these two unresolved lines are well measured. 
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